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and compare it with the corresponding Landau gauge propagator, showing that for both the rele- 
vant IR mass scale coincides. We also report preliminary results on Coulomb gauge ghost form 
factor and quark propagators and give a comment on the gluon propagator's strong coupling limit. 
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Introduction 

Yang-Mills theories wave functional studies in Coulomb gauge have enjoyed a renewed inter- 
est in the last years [1, 2, 3, 4, 5], since they bypass an explicit construction of the physical Hilbert 
space [6]. The Gribov-Zwanziger confinement scenario plays a key role in this context [7, 8] and 
its predictions for the static two pint functions need therefore to be explicitly verified. We re- 
port here on our lattice calculation program regarding static gluon and ghost, the relation between 
Coulomb and Landau gauge and the gluon strong coupling limit. We also show preliminary results 
for the quark propagator in Coulomb gauge, which should relate the Gribov-Zwanziger mechanism 
to chiral symmetry breaking. 



1. Gluon propagator 

As shown in [9, 10], after taking care of subtle renormalization issues related to the energy 
dependence, the SU (2) static gluon propagator in Coulomb gauge D(\p\) agrees within numerical 
precision with Gribov's formula [7], with an IR mass M = 0.856(8)GeV and vanishing anomalous 
dimension in the UV. In [11] the relation between Coulomb and Landau gauge gluon propagators 
was futher investigated. It was shown that both in 3 and 4 dimensions the two propagators coincide, 
up to a simple rescaling p(p) of the momentum variable, shown in Fig. 1. The rescaling function 
in 3 dimensions is a simple interpolation between two constants, since both propagators have no 
anomalous dimension. In 4 dimensions only the Landau propagator has an anomalous dimension; 
the function p behaves accordingly. When redrawn as a function of pp(p), as shown in Fig. 2, the 
Landau gluon propagator Di{p) coincides with the equivalent Coulomb one, Dc(p) = \p\~ l D(\p\); 
both show a massive IR behavior. Whether dimension 2 condensates [12, 13] could be related to 
such IR mass is still an open question. 
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Figure 1: Momentum rescaling function between Coulomb and Landau gauge propagators. 
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Figure 2: Landau and Coulomb gluon propagator vs. equivalent momenta in 3 (left) end 4 dim. (right). 

2. Strong coupling limit 

The massive IR behavior in Landau gauge, found in all standard lattice calculations [14, 15], 
is in disagreement with the Kugo-Ojima picture [16]. Some authors [17, 18] have supposed that 
discretization and/or Gribov copies effects might be responsible for the discrepancies; strong cou- 
pling calculations have been put forward as evidence. However going away from the continuum 
limit one goes through bulk transitions [19] and unphysical effects should be expected, so that pre- 
dictions arising from (5=0 results should be taken with extreme care. As an example, we show 
the temporal gluon propagator in Coulomb gauge Dq(\p\), which should be equivalent to the static 
potential [8]. Even without going to anisotropic lattices [20], which take care of scaling violations 
[9, 10, 21], its diverging behavior, as expected for a confining theory, has been long established 
in the continuum limit. In Fig. 3 Do(|p|) is given in lattice units; the physical behavior seems to 
be completely lost at strong coupling. The only way out would be to accept that the lattice scale 
diverges in the j8 — > limit. All data would then collapse to a single point at (0,°°) when expressed 
in physical units. The physical information contained in such result is of course minimal. 
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Figure 3: Coulomb gauge temporal gluon propagator Dq(\p\) in the strong coupling limit. 
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3. Ghost form factor 

Functional methods [3, 4] provide sum rules which must be satisfied by the static gluon 
propagator and ghost form factor IR exponents K g i + 2K g h = 1 and UV anomalous dimensions, 
Y g i + 2j g i, = 1, and give precise predictions for their value, K g i = — 1, K g h = 1, J g i = and y g h = 1 /2. 
For the static gluon propagator these have been already verified [9, 10, 11]. Extending previous 
SU(2) calculations [22], new preliminary results for the ghost form factor on a 32 4 lattice, 

shown in Fig. 4, seem to confirm such prediction. Indeed, we find that the data are very well fitted 
in the whole momentum range by: 



d(\p\)ocJ t ^+]og-\e + ^), (3.1) 
y p L m z 

with m = 0.32(2)GeV and # 2 /d.o.f. = 0.8. Calculations on larger lattices are however necessary 
and under way to confirm this result [23]. 
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Figure 4: Ghost form factor rf(|p|) on a 32 4 lattice. 



4. Quark propagator 

Although not directly related to the Gribov-Zwanziger confinement scenario, the quark propa- 
gator in Coulomb gauge provides important physical information on chiral symmetry breaking and 
dynamical mass generation. We report here preliminary results obtained in SU (3) for quenched 
12 3 x 24 lattices, using asqtad staggered fermions with a bare mass of mi, ~ 212 MeV. Results 
using unquenched configurations and masses closer to the chiral limit will be published soon [24]. 
We have followed [25, 26, 27], adapting their method and notation from Landau to Coulomb gauge, 
and explicitly inverted the Dirac operator on each configuartion. Our first non-trivial results is that 
the coefficients of the mixed term fa p, which are in principle allowed in Coulomb gauge beyond 
tree level, vanish within numerical precision. We find then that the most general form for the 
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propagator and its inverse reads: 



S(p,po) = -ifo^t(p,po)-ip^(p,po) + & m (p,Po) 
S~ l (p,p Q ) = ip'oA t (p,po) + ip A s (p,p ) + B m (p,p ) , 



where A t , A s and B m can be easily re-expressed in terms of g/ t , srf s and 
obtained from the simulations. If renormalizable, one should expect: 

S(p,Po) = Z(p,p ) 
1 



(4.1) 

, which are explicitly 



S (p,Pa) 



if + ipo a(p,po)+M(p,po) 
[ip + ifioa(P,Po) +M(p,p )} 



(4.2) 



Z(p,Po) 

and Z(p,po), a(p,po) and M(p,po) can be again expressed in terms of stf t , s/ s and 3§ m \ Z(p,po) 
should be multiplicatively renormalizable while the functions a(p,po) and M(p,po) should be 
cut-off independent. Fig. 5 shows the former for two different values of the coupling; it is clear 
that the full quark propagator S(p,po) is not renormalizable. This seems to be a general pattern: 
just like for the full Coulomb gauge gluon propagator [9, 10, 11], the po dependence needs to be 
treated separately. Turning to the static quark propagator S~ l (p) = Jd poS~ l (p,po), being a an 
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Figure 5: Temporal component a(k,ko) of the quark propagator in Coulomb gauge, 
even function of po, the temporal component cancels for parity and we obtain: 

S(P) = Z{p) 7 



where: 



M(p) 

m 



s-\p) 



ip + M{p) 
^ [if+M{p)\ , 



m 

Jdp B m (p,po) Jdpo^ m (p,p ) 



Jdp A s (p,p ) fdp g/ s (p,p ) 

r , ,V S = [P 2 + M2 (P)] I dpo^ s (p,Po). 
Jdp A s (p,po) A J 



(4.3) 



(4.4) 
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Figure 6: Static quark self energy Z(k) and running mass M{k). 



Z(p) is proportional to the quark dispersion relation and should be multiplicative renormalizable, 
while the running mass M(p) should be renormalization group invariant. The two functions are 
shown in Fig. 6, where Z(p) has been rescaled for different /3 while the values for M(p) have 
been directly taken from Eq. 4.4 without modifications; the light discrepancy is simply due to the 
different renormalized current mass for different cutoffs given a fixed bare mass m^. A cone cut 
has been applied to both data sets. 



Conclusions and Outlook 

We have reported on the progress in our program to determine two point functions in Coulomb 
gauge. The static gluon propagator and ghost form factor confirm down to the available IR mo- 
menta the Gribov-Zwanziger confinement scenario and the predictions from the hamiltonian func- 
tional approach. We also show that, as for the gluon, only the static component of the quark 
propagator is renormalizable. Extension to higher volumes and unquenched configurations are cur- 
rently under investigation [23, 24]. We also comment on strong coupling limit results for the gluon 
propagator. We show that, as expected, physical interpretations in this limit must be taken with 
extreme care. 

Acknowledgments 

We thank D. Campagnari and P. Watson for helpful discussions. This work was partly sup- 
ported by DFG under the contract Re856/6-3. 

References 

[1] A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D65, 025012 (2002), [hep-ph/0107078]. 

[2] C. Feuchter and H. Reinhardt, Phys. Rev. D70, 105021 (2004), [hep-th/0408236]. 

[3] W. Schleifenbaum, M. Leder and H. Reinhardt, Phys. Rev. D73, 125019 (2006), 
[hep-th/0605115]. 



6 



Propagators in lattice Coulomb gauge and confinement mechanisms 



Giuseppe Burgio 



[4] D. Epple, H. Reinhardt and W. Schleifenbaum, Phys. Rev. D75, 04501 1 (2007), 
[hep-th/0612241]. 

[5] M. Quandt, H. Reinhardt and G. Burgio, Phys. Rev. D81, 065016 (2010), [arXiv: 1001 . 3699]. 

[6] G. Burgio, R. De Pietri, H. A. Morales-Tecotl, L. F. Urrutia and J. D. Vergara, Nucl. Phys. B566, 
547-561 (2000), [hep-lat/9906036]. 

[7] V. N.Gribov, Nucl. Phys B139, 1 (1978). 

[8] D. Zwanziger, Nucl. Phys. B485, 185 (1997), [hep-th/9603203]. 

[9] G. Burgio, M. Quandt and H. Reinhardt, Phys. Rev. Lett. 102, 032002 (2009), 
[arXiv: 0807 . 3291]. 

[10] G. Burgio, M. Quandt and H. Reinhardt, PoS CONFINEMENTS, 051 (2008), 
[arXiv: 0812 . 3786]. 

[11] G. Burgio, M. Quandt and H. Reinhardt, Phys. Rev. D81, 074502 (2010), [arXiv : 9 1 1 . 5 1 1]. 

[12] G. Burgio, F. Di Renzo, G. Marchesini and E. Onofri, Phys. Lett. B422, 219 (1998), hep-ph/9706209. 

[13] P. Boucaud et al., JHEP 04, 006 (2000), hep-ph/0003020. 

[14] I. L. Bogolubsky et al., Phys. Rev. D74, 034503 (2006), [hep-lat/0511056]. 

[15] I. L. Bogolubsky et al., Phys. Rev. D77, 014504 (2008), [arXiv: 0707 .3611]. 

[16] T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66, 1 (1979). 

[17] A. Sternbeck and L. von Smekal, Eur. Phys. J. C68, 487 (2010) [arXiv : 8 1 1 . 4 3 0]. 

[18] A. Maas, J. M. Pawlowski, D. Spielmann, A. Sternbeck and L. von Smekal, Eur. Phys. J. C68, 183 
(2010) [arXiv: 0912 . 4203]. 

[19] G. Burgio, PoS LAT2007, 292 (2007), [arXiv: 0710 .047 6]. 

[20] G. Burgio et al. (TrinLat), Phys. Rev. D67, 114502 (2003), [hep-lat/0303005]. 

[21] Y. Nakagawa, A. Nakamura, T. Saito and H. Toki PoS LAT2009, (2009), [arXiv: 911 . 2550]. 

[22] M. Quandt, G. Burgio, S. Chimchinda and H. Reinhardt, PoS LAT2007, 325 (2007), 
[arXiv: 0710 . 0549]. 

[23] M. Quandt, G. Burgio and H. Reinhardt, In preparation, (2010). 

[24] M. Schrock, G. Burgio, M. Quandt, and H. Reinhardt, In preparation, (2010). 

[25] J. I. Skullerud and A. G. Williams, Phys. Rev. D63, 054508 (2001), [hep-lat/0007028]. 

[26] P. O. Bowman, U. M. Heller and A. G. Williams, Phys. Rev. D66 014505 (2002), 
[hep-lat/0203001]. 

[27] P. O. Bowman et al. Phys. Rev. D71 054507 (2005), [hep-lat/0501019]. 



7 



